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Abstract 

We construct exact vortex solutions in 3 + 1 dimensions to a theory which is 
an extension, due to Gies, of the Skyrme-Faddeev model, and that is believed to 
describe some aspects of the low energy limit of the pure SU{2) Yang-Mills theory. 
Despite the efforts in the last decades those are the first exact analytical solutions 
to be constructed for such type of theory. The exact vortices appear in a very parti- 
cular sector of the theory characterized by special values of the coupling constants, 
and by a constraint that leads to an infinite number of conserved charges. The 
theory is scale invariant in that sector, and the solutions satisfy Bogomolny type 
equations. The energy of the static vortex is proportional to its topological charge, 
and waves can travel with the speed of light along them, adding to the energy a 
term proportional to a [/ (1) Noether charge they create. We believe such vortices 
may play a role in the strong coupling regime of the pure SU{2) Yang-Mills theory. 
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The Skyrme-Faddeev (SF) model [T] has been proposed long ago as a four dimensional 
field theory that can present finite energy knot solitons carrying a topological charge given 
by the Hopf map 5*^. Despite many efforts in the last decades no analytical non- 

trivial solutions have been constructed. The confirmation of that expectation however 
came some years ago with the construction of numerical solutions employing powerful 
computer resources [2], El The interest in the SF model has increased since then 
and many applications have been proposed including Bose-Einstein condensates |S] and 
superconductors [6J. In addition, it has been conjectured [7] that the SF model describes 
the low energy limit of the pure (no matter) SU{2) Yang-Mills theory, with the knot 
solitons being interpreted perhaps as glueballs, or some non-trivial vacuum configurations. 
The conjecture is based on the Cho-Faddeev-Niemi decomposition [8], [7] of the SU{2) 
Yang-Mills field in terms of an abelian gauge field, a triplet n of scalars living on S*^ 
(n^ = 1), and two other degrees of freedom which can be expressed by two real scalar 
fields. There has been several controversies about the validity of that conjecture P, HO] 
but a calculation carried out by Gies [H] has indicated that it may indeed be the case 
provided additional quartic terms are added to the effective Lagrangian. 

The purpose of this paper is to present the first exact analytical solutions in 3 + 1 
dimensions for the extension of the Skyrme-Faddeev model defined by the Lagrangian 

£ = d^n ■d^'n-^ {d^n A d^nf + - (S^n ■ d^nf (1) 

where n is a triplet of real scalar fields taking values on the sphere S*^, M is a coupling 
constant with dimension of mass, and 13 are dimensionless coupling constants. We show 
that ([T]) possesses a Bogomolny type sector admitting an infinite number of conserved 
currents, and when the coupling constants satisfy (3 = 1, there exists an infinite class 
of exact solutions. Among those there are static vortex solutions as well as vortices with 
waves travelling along them, and enhancing their stability. 

The first two terms of ([1]) correspond to the original SF model [H H], and the third 
one agrees with the term found by Gies in his calculations [TTJ. A similar extension of 
the SF model has been considered in [12], but with some important differences. First the 
model was defined in three space dimensions only and the signs of the second and third 
terms in ([1]) were chosen to keep the static Hamiltonian positive definite. The existence 
of the exact solutions we present in this paper depends crucially on the fact that and /3 
have the same sign, and that agrees with the signature of the one-loop Wilsonian effective 
action for the SU{2) Yang-Mills theory obtained in [llj. Despite of that signature, our 
solutions belong to a sector of the theory where the energy is positive definite. Solutions 
for the conformally invariant version of the model, namely the theory ([T]) without the first 
term (M = 0), have been constructed and will be given in a separate paper jl3] . 
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The solutions we construct in this paper can be best explained by using the stereo- 
graphic projection of 5*^ in terms of a complex scalar field u, i.e. 



n 



One then obtains that 



(u + u*, -i {u-u*),\u p -l) / {l+\u 



n ■ {d^n A dt,n) - 

{d^n ■ d^n) -- 
Therefore, the Lagrangian ([T]) becomes 



2i 



(1+ I U 

d„u d^u* 



{i+\u\^y 



1+ I u P)2 + e2 



^2 _ ^^ {d,u d^u* 



\2n 



and the Euler-Lagrange equations following from (jl]), or ([1]), reads 

(l+ I u |2) d^'lC^ - 2 u* /C^ d^'u = 
together with its complex conjugate, and wher^ 



4 [(1 - /5e2) (d^ud^u*) d^u - {d^ud^'u) d^u*] 



[1+ I u 



(2) 



(3) 



(4) 



(5) 



(6) 



It was shown in using a generalization of the zero curvature condition to higher 
dimensions, that many field theories admit an infinite number of conserved charges when 
some special constraints are imposed. In theories with target space S"^ like ([T]) the relevant 
constraint is given by 

(7) 



d^u d^u = 



When ([7]) is imposed the theory ([T]) possesses the infinite set of conserved currents given 
by 



6G 



6G 



where G is any function of u and u*, but not of their derivatives, and /C^ is obtained from 
(E]) by imposing (I7j), i.e. 



^ ^^62 (l+\u 



(9) 



^Notice that this theory has a set of trivial plane wave solutions given by m = Ae f ^ , with the 
wave vector satisfying either fc^ = or fc^ = — -^^ - 



^^^i' Plane wave solutions for the SF model were 



considered in [T4l. 
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The conservation of ([8]) follows from the equations of motion, which now read d^lC^^ = 0, 
and the two identities IC^^d^u = and Im = 0. In the context of SF model 

those currents were first constructed in [T7] . 

Notice that the constraint ([7]), in Minkowiski space-time, can be cast into the form 

[(^1 + I 82) U] -td2)u] = - [(^3 + 9o) U] [(^3 - 9o) U] (10) 

If one now chooses the coupling constants in ([T]) to satisfy 

(3e^ = l (11) 

one gets that eqs. of motion ([5]) reduce to d^u = 0, or 

(^1 + 1 82) (9i - z ^2) n = - (^3 + ^o) (^3 - do) u (12) 

In a Euclidean space-time the same relations hold by replacing the time coordinate x° by 
an imaginary time —ix^, and so ^ id^. 

Of course the equations (fTOj) and (fT2|) are solved by field configurations satisfying 

{di + 2 £1 ^2) M = and (^3 + £2 f^o) w = (13) 

where = ±1, are signs chosen independently. In addition, equations ( fT3l) are satisfied 
by configurations of the formal 

u = V {z) w {y) (14) 

where z = + isix^ and y = x^ — 62 x'^, with v {z) and w {y) being arbitrary regular 
functions of their arguments. Notice that if u satisfies f|T3l) . so does any regular functional 
of it, J-" (u). Indeed, by taking JF to be the logarithm one observes that the ansatz ( fT4l) is 
mapped into u = v (z) + w (y). See [19] for similar discussions in 2 + 1 dimensions. 

Notice that by imposing the constraint (171) and the condition ( iTTl) the Lagrangian (jlj) 
and equations of motion ([5]) becomes those of the CP^ model, i.e. 

Ccpi = 4 (l+\u\^) d\-2u* {d^'uf = (15) 
(1+ I M P) 

Therefore, the solutions of (ITO!) and (fT2|) also solve the CP^ model in 3 + 1 dimensions. 
Therefore, the Skyrme-Faddeev-Gies model ([T]) shares a common submodel with CP^, 
defined by the equations {d^uf = 0, and d'^u = 0, together with the condition (ITT]) . In 
fact, the first equation in (fT3l) correspond to the Bogomolny equation, or equivalently the 



■^We are grateful to Prof. D. Fairlie for pointing to us that a large class of implicit solutions for the 
wave equation and the constraint ([7]) is given in [18], and that they are related to the solutions ^T^ ■ 
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Cauchy-Riemann equations, leading to static solutions of that model in 2 + 1 dimensions 
[20] . The relation between the Bogomolny equation and the constraint ([71) was already 
pointed out in J6] when the generalization of the zero curvature condition was applied 
to the CP^ model in 2 + 1 dimensions. 

The structure of the solutions is now clear. The function v (z) corresponds to the static 
lumps of the CP^ model in 2 + 1 dimensions. The function w (y) corresponds to waves 
traveling along the direction with the speed of light. In a Euclidean space-time w (y) 
corresponds instead to a second copy of the CP^ lumps. If one wants static solutions in 
Minkowiski space-time, then it has also to be independent. Therefore, the CP^ lump 
gets a string like shape in three space dimensions and it corresponds in fact to a static 
vortex solution for the Skyrme-Faddeev-Gies model (alternatively for the CP^ too, in 
three space dimensions). One can take v (z) to be any finite energy lump solution of the 
two dimensional CP^ model. For instance, one can take f ~ z", and w (y) = 1 in f|T4l) . 
Using polar coordinates on the x^ plane, i.e. x^ +ieix'^ = p e*^^ we obtain the static 
vortex 



where n is an integer, and a an arbitrary parameter with dimension of length. Config- 
urations with several vortices all parallel to the x^-axis can be obtained from the multi 
lumps of the CP^ model. Numerical and approximate vortex solutions for the SF model 
have been considered in [Ti] . 

One can dress such vortices with waves traveling along the x^ axis. There are several 
ways of doing it. One way of keeping the energy per unit of length finite is to take w (y) 
in (JT4l) of the plane wave form, leading to the vortex 



where k is an arbitrary parameter with dimension of (length) ^ . Vortices with waves have 
been found in a different context in P^ . 

If one evaluates the Hamiltonian density for the theory ([T]) and then imposes the 
conditions (ITIl and ffTTl) one obtains 



where the subindex c means we are imposing ffTOl) and ffTTl) . Notice that it coincides with 
the Hamiltonian density for the CP^ model. We can write that as 




(16) 




(17) 




(18) 



He 



diu + ieid2U |^ + | d^u + e2dou 



(1+ I U |2) 
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+ iei {diud2U* — d2udiu*) — £2 {d^udou* + doudsu*)] 
{diud2U* - d2udiu*) 



> AM' 



iSr 



e2 



(1+ |« |2) 

{d^udou* + doud^u*) 
(1+ U P)' 



(19) 



The bound on the energy density is clearly saturated by the solutions of f|T3|) . In fact, for 
those solutions one gets 



Hc = 8M^ (1 d,u p + I dyu |2) / (1+ I M 



(20) 



The same analysis can be done in the Euclidean case, with the same results, replacing the 



-I X 



Hamiltonian density by the Lagrangian density, and - 

The energy density for the vortex solutions (1161) and (IT7|) is then given bjfl 



vortex 



8M' 



2n 



(21) 



Therefore, by integrating on the x^x"^ plane we get the energy per unity lenght. For the 
static vortex (fT6l) it is 



Stat. vortex 



n 



(22) 



The integer n is the topological charge associated to the vortex, and defined as the winding 
number of the map from any circle on the x^x^ plane, centered at the x^-axis, to the circle 
m/ I M I on target space. So, the energy per unit length is proportional to the topological 
charge, a characteristic of solutions saturating the Bogomolny bound. 

The energy per unit length for the time-dependent vortex f[T7|) diverges for n = ±1, 
and for | n |> 1 it is given by 



vortex / wave 



SttM^ 



n I +k'a'I{\ n |) 



(23) 



where / {n) = :^ E (~") T (~"); with E being the Euler's Gamma function. Notice that 
I [n) is a monotonically decreasing function of n and / (n) l/n, as n ^ 00. Therefore, 
the plane wave contribution to the energy per unit length decreases as | n | increases. 

Notice that a necessary condition for the charge densities Jq associated to ([8]) to be 
non- vanishing is that the field u should be time dependent. Therefore, all the charges 



Notice that such energy density depends on | n |, and not on the sign of n, because 



4 



(i+(p/a) 
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associated to ([8]) vanish when evaluated on the static vortex solutions (fT6|) . However, 
the vortices (|T71) . with waves travelling along them, are time dependent solutions and so 
can have non- vanishing charges associated to ([8]). Consequently that infinite number of 
conserved charges introduce selection rules which protect them against decay into lower 
energy vortex solutions, and so give them some stability. However the question of stability 
of such solutions under small perturbations is a much more complex issue, and deserves 
further studies. We should point out however that the static Hamiltonian associated to 
([1]) is positive definite for > 0, < 0, /? < and [3 > 1. Therefore, our static 
vortices, saturating the Bogomolny bound (see (fT9l) ). have the smallest possible energy 
in such theory, for that range of coupling constants. That fact makes us believe that for 
sufficiently large M^, the vortices should be stable under small perturbations. 

If one chooses in ([8]), the functional as G = — ^4 (1+ | u p) ^, and if one takes into 
account (ITT!) one gets that the corresponding conserved current is 



which corresponds to the Noether current for the theory ([T]) associated to the symmetry 
u u, after the conditions ([7]) and (ITT!) are imposed. Evaluating the charge per unit 

length for the solution (ITTI) one gets 



Therefore, the contribution to the energy per unit length fl23|) coming from the plane wave 
is proportional to the U{1) Noether charge it gives to the solution. 

The spectrum of energy we have got is therefore very similar to that of the so-called 
Q-lumps constructed in [2T] for the 2 + 1 dimensional CP^ model with a potential. See 
[22] for more details, [23] for a related model, and also [24J for Q-balls solutions on S^. 
However, the mechanisms involved here are quite different. We have a four dimensional 
theory, without a potential term depending on | m | only, and our vortices are not spinning 
but instead have waves traveling along them with the speed of light. In addition, we again 
point out that our results are true for the four dimensional CP^ model without a potential, 
but subjected to the constraint ([7]). 

The theory ([1]) is invariant under the internal 5*0(3) symmetry of rotations of the n 
field. It then possesses three conserved Noether currents associated to such symmetry. 
One of them is that associated to the symmetry u e^^u, discussed above. However, 
when the contraint ([7]) is imposed the reduced model gets an infinite number of conserved 
currents given by ([H]). Such currents can not be of the Noether typ^ because the reduced 

'^Except for those three currents associated to the global 50(3) symmetry of ([T]), contained in ([8]). 



J^ = -iA M- 



2 [u d^u* - u* d^u] 



(24) 




(25) 



6 



system of equations, namely ff^u = and (S^m) = 0, does not possess a Lagrangian 
and so the Noether theorem does not apply. Such currents can be understood as hidden 
symmetries of the reduced model by constructing them using a zero curvature represen- 
tation as explained in section 6.1.1 of |16]. They are associated to an infinite dimensional 
non semi-simple Lie algebra made of the semi-direct product of SU{2) with abelian ide- 
als transforming under infinite dimensional representations of SU{2). It is important to 
understand the role such symmetries may have in the extension of the Skyrme-Faddeev 
model given by the theory ([1]). 

Another point concerning the enlargement of symmetries is the following. Even though 
the theory ([1]) is not scale invariant, the conditions ([7]) and (11 II) lead to a sub- model which 
equations of motion are scale invariant. The condition (du)"^ = is conformally invariant 
but d'^u = is only scale invariant in four dimensions. 

Finally, we would like to comment on some important points of our construction which 
relate to the pure SU{2) Yang-Mills theory. Notice that by using the Cho-Faddeev-Niemi 
decomposition [8], [7] of the SU{2) gauge field A^, and assuming that at low energies all 
degrees of freedom get frozen except for those associated to the n field, then = d^n/\n. 
Our vortices are solutions of the constraints d^ud^u = 0, and also of d'^u = 0. Those 
equations imply that the field n satisfies d'^fl + ((9n)^ n = 0. Therefore the connection, 
evaluated on the vortices, satisfies the gauge fixing condition 

d^A^ = (26) 

Another point is that, without the use of any condition (equations of motion or constraint), 
and working in a Euclidean space-time, it follows that (see eq. ([3])) 

1^ . 1^ = Hc/M^ (27) 

where He is given in (fTSl) . Therefore as explained in (fT9l) . the equations (fT3l) then imply 
that our vortex solutions minimize locally the quantity A^- A^. In fact, its integral on the 
plane x"^ has the smallest possible value when evaluated on the vortices. Such quantity 
is used as a gauge fixing condition and it is believed to signal non trivial structures in 
gauge theories [26] . 

In order to have the condition (ITT]) satisfied we need of course the coupling constants 
13 and to have the same sign. That implies that the second and third terms in ([T]) 
have opposite signs. We are then in a situation which differs from that considered in [T2] 
where the signs of the three terms in ([1]) were chosen to give positive contributions to the 
static energy. See also |25j for a discussion on solutions of that type of model. We are 
however in agreement with the situation in [H]. Indeed, consider the action ([T]) in four 
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dimensional Euclidean space and take (3 and to be negative. We then have the same 
sign structure of the effective action given in eq. (14) of [H]. In fact, the condition ( ITIi) 
corresponds to the following relation among the parameters of [11] 



where is the gauge coupling constant, a a gauge fixing parameter, A and k are the UV 
and IR cutoffs respectively. The values for a in [11] are such that (4 — a) /a is positive, 
and so (!28l) is indeed consistent with k < A. However, the relation (!28|) has to be read 
in its proper context, i.e. the physical limit corresponds to A; — 0, and the perturbative 
results of [TT] may not be applicable in that regime. It would be very interesting to in- 
vestigate if the condition (fTTl) remains consistent with the renormalization group flows in 
the IR limit. That could clarify if the solutions calculated in this paper can play a role 
in the low energy limit of the pure SU (2) Yang- Mills theory. 
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